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transionospheric  propagation  channel.  It  is  shown  that  aperture  averaging 
can  reduce  observed  signal  power,  increase  observed  decorrelation  time  and 
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anisotropic  electron  density  irregularities  that  are  elongated  in  the 
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field  line  with  respect  to  the  direction  of  propagation  is  arbitrary.  This 
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SECTION  1 
INTRODUCTION 


Large  high-gain  antennas  are  used  in  many  radar  and  communica¬ 
tions  applications  to  increase  the  energy  collected,  to  increase  angular 
accuracy,  and  to  provide  protection  against  jamming.  If  the  wavefront  at 
the  antenna  aperture  experiences  scintillation,  large  apertures  can  act  to 
average  the  signal  and  thereby  modify  the  observed  signal  properties. 
Aperture  averaging  is  important  for  large  lenses  at  optical  wavelengths 
(Tatarskii,  1971)  but  is  not  generally  important  at  satellite  frequencies 
under  ambient  or  naturally  perturbed  ionospheric  conditions. 

However,  communications  and  radar  systems  from  VHF  to  SHF  that 
utilize  trans-ionospheric  propagation  geometries  can  experience  severe 
Rayleigh  fading  after  a  high  altitude  nuclear  detonation.  This  increased 
level  of  scintillation  is  caused  by  the  creation  of  great  quantities  of 
ionization  which  evolves  into  large  irregular  structure  that  takes  the 
form  of  striations  or  filaments  aligned  along  the  earth's  magnetic  field. 
Since  the  velocity  of  propagation  is  dependent  upon  the  local  index-of- 
refraction,  different  portions  of  a  wavefront  propagate  at  different 
velocities  in  a  striated  region.  Thus  a  once  plane  wavefront  can  become 
distorted  with  different  portions  traveling  in  slightly  different  direc¬ 
tions.  As  the  wave  propagates  farther,  diffraction  or  angular  scattering 
causes  constructive  and  destructive  interference  which  introduces  fluctua¬ 
tions  in  amplitude,  phase,  and  angle-of-arrival. 


For  cases  of  severe  scintillation  where  the  signal  varies  across 
the  area  of  the  receivinq  aperture,  the  effect  of  the  aperture  can  be 
significant.  It  is  well  known  that  an  aperture  antenna  acts  to  cut  off 
enerqy  that  is  incident  at  off-boresight  anqles  where  the  antenna  gain 
function  is  reduced.  This  effect  may  also  be  viewed  as  the  result  of 
averaging  or  coherent  processing  of  the  electromagnetic  field  received  at 
the  aperture  location.  In  this  report  the  effects  of  aperture  averaging 
are  analytically  calculated  for  the  situation  where  a  Gaussian  antenna 
beam  is  used  by  the  receiver.  It  is  shown  that  aperture  averaging  has  a 
greater  effect  on  antenna  measurements  of  mean  time  delay  and  time  delay 
jitter  than  on  measurements  of  decorrelat ion  time  or  scattering  loss. 
Simple  analytic  expressions  are  given  for  all  these  quantities  in  terms  of 
the  geometry  of  the  propagation  path  and  the  severity  and  structure  of  the 
ionization  irregularities. 

To  obtain  results  applicable  to  a  qeneral  geometry  of  the  line- 
of-sight  relative  to  the  field  aligned  ionization  structure,  it  is  neces¬ 
sary  to  obtain  an  analytic  solution  for  the  two-position,  two-frequency 
mutual  coherence  function  for  spherical  wave  propagation  throuqh  a  thick 
layer  of  anisotropic  electron  density  irregularities.  It  is  assumed  here 
that  strong  scattering  conditions  prevail  and  that  the  quadratic  approxi¬ 
mation  to  the  phase  structure-function  is  valid. 

This  approximation  was  used  by  Sreenivasiah  et  al.,  (1976)  and 
by  Sreenivasiah  and  Ishimaru  (1979)  for  the  cases  of  plane  wave  and  beam 
wave  propagation  in  homogeneous  turbulence.  More  recently  the  two-posi- 
tion,  two-frequency  mutual  coherence  function  was  obtained  for  spherical 
wave  propagation  using  the  extended  Huygens-Fresnel  principle  (Fante, 
1981).  Although  the  quadratic  structure-function  approximation  can  some¬ 
times  lead  to  difficulties  (Wandzura,  1980)  it  is  appropriate  for  the 
two-frequency  mutual  coherence  function  but  not  for  calculation  of  higher 
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moments  of  the  field  (Fante,  1980).  Fante  (1981)  discusses  the  accuracy 
of  the  quadratic  structure-function  approximation  for  the  case  of  atmo¬ 
spheric  turbulence  with  a  Kolmogorov  power  spectrum.  He  has  found  that 
the  accuracy  is  a  function  of  the  irregularity  power  spectrum  and  of  the 
strength  of  the  turbulence  (Private  Communication,  1982),  with  accuracy 
increasing  for  stronger  scattering. 

This  work  is  a  genera  1  izat ion  of  an  earlier  calculation  (Kne,op, 
1982)  valid  only  for  isotropic  turbulence.  As  in  the  former  study,  the 
results  here  are  specialized  to  the  case  of  a  thin  phase-screen  approxima 
tion  to  the  thick  scattering  medium.  With  this  simplification  analytic 
results  are  obtained  for  the  received  impulse  response  function  to  a 
transmitted  power  delta  function.  Then  results  may  easily  be  determined 
for  the  mean  time  delay  and  time  delay  jitter  for  stronq,  anisotropic  tur 
bulence  in  the  thin  phase-screen  approximation.  Results  for  these  quant i 
ties  in  the  thin  phase-screen  approximat ion  have  previously  been  shown 
(Knepp,  1982)  to  closely  approximate  the  results  for  a  thick  scattering 
layer  for  the  case  of  isotropic  turbulence. 

The  effects  of  a  Gaussian  receiving  antenna  on  measurements  of' 
received  power,  decorrelation  time  (or  distance),  mean  time  delay,  time 
delay  jitter  and  coherence  bandwidth  are  determined.  It  is  shown  that 
aperture  averaging  can  reduce  observed  signal  power,  increase  observed 
decorrelation  time  and  be  a  significant  factor  in  reduced  observed  time- 
of-arrival  jitter  at  the  antenna  output. 


SECTION  2 
FORMULATION 


Consider  a  monochromatic  spherical  wave  E("p,z,u),t)  which  origi¬ 
nates  from  a  transmitter  located  at  (0,0,-zt)  and  propagates  in  free- 
space  in  the  positive  z  direction  where  it  is  incident  on  an  ionization 
irregularity  layer  which  extends  from  0  z  £  L.  After  emerqing  from  the 
layer  at  z  =  L,  the  wave  then  propagates  in  free-space  to  a  receiver 
located  at  (0,0,zr).  This  thick  layer  geometry  is  shown  in  Figure  1. 

As  the  wave  propagates,  its  phase  substantially  behaves  as  (-i<k>z+iwt)  so 
write 


E(p,z,w,t)  =  U(p,z,w)  exp{i(u)t-/<k(z')>dz') } 


(1) 


Figure  1.  Propagation  of  signals  through  a  disturbed 
trans ionospheric  channel. 
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where  <k(z)>  is  the  mean  wave  number  given  by 

0)  V2  ,  ,  l/?- 

<k(z)>  =  “  (1  -  Ne/nc)  =  k  (1  -  k2/k2)  (2) 

where  c  is  the  speed  of  light  in  a  vacuum,  Ne  is  the  mean  ionization 
density,  nc  is  the  critical  electron  density  and  is  related  to  the  clas¬ 
sical  electron  radius  r^  by  nc  =  ir/(X2re),  (re  =  2.8?*10-15  m). 

It  has  been  shown  that  U(T>,z,w)  satisfies  the  parabolic  wave 
equation  under  the  Markov  approximation,  (Tatarskii,  1971:  Yeh  and  Liu, 
1977) 

AN  (o2 
_ e  _p 

z  an  ,  <N  >  uj2 

Vi  IJ  -  2i<k(z)>  2L  -  <k2(z)>  _J? _  U  =  0  (3) 

“  3Z  (l-Wp/u)2) 

where  ANe  is  a  small  variation  in  the  ionization  level  and  wp  =  kpc 
is  the  circular  plasma  frequency  of  the  background  ionization.  The 
exponential  time  dependence  has  been  suppressed  and  it  has  been  assumed 
that  up  <  w,  otherwise  signal  attenuation  would  be  the  dominant  effect. 

2.1  MUTUAL  COHERENCE  FUNCTION 

Now  in  the  case  that  the  transmitted  waveform  is  no  lonqer  a 
monochromatic  wave,  but  can  be  expressed  as  a  waveform  modulated  on  a 
carrier,  the  two-position,  two-frequency  mutual  coherence  function  r  is  of 
interest,  r  is  important  for  the  calculation  of  pulse  propagation  in  a 
random  medium  and  it  serves  as  a  basis  from  which  to  calculate  the 
important  power  impulse  response  function  and  its  moments.  Under  the 
Markov  approximation,  r  satisfies  the  following  equation  (Tatarskii,  1971; 
Yeh  and  Liu,  1977). 
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ar 

_ + 

3z 


1 


2kik2 


(k2?2  -kiV2  )r 

■M  J-2 


•  -g  [2k iBik2B2A( P1-P2) 


(k12e12+k22B22)A(0) ]r  =  0  (4) 


where 


r  -  <u(xi,yi,z,(i)i)U*(x2,y2,z,w2)> 


and  ki  and  k2  are  the  wavenumbers  at  two  frequencies  and  o>2, 
respectively,  and  V2  is  the  two-dimensional  Laplacian 

j-i 


v2 

M 


+ 


3 


2 


with  a  similar  definition  for  v2  .  Following  Veh  and  Liu  (1977)  let 


w  2/ U) , 2 

P  1 

1-01  2/oi,2 
p  1 


(5) 


with  a  similar  definition  for  02. 


The  function  A(p)  is  the  integral  of  the  autocorrelation  func 

tion  of  electron  density  fluctuations,  B$,  in  the  direction  of  propa- 
qat ion 


OO 

A(pi-Pz)  =  /  Cp i -"p2 , z  ' ) dz ' 


(6) 
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where 


Pi"P2  =  ( x i~x  2.y i-y 2) 


£  »  _^e 

<NP> 


so  that 


_  _  "  ’^.*(pi"P2) 

A(pi-p2)  =  2 »  /  /  e  *c(K^,Kz=0)d2KJ 


(7) 


where  is  the  power  spectrum  of  electron  density  fluctuations.  Equa¬ 
tions  6  and  7  depend  upon  the  validity  of  the  Markov  approximation  where 
it  is  assumed  that  the  electron  density  fluctuations  are  delta-correlated 
in  the  direction  of  propagation  (Fante,  1975).  That  is 


^ Pl*P2,z i'z2)  =  A( pj-p2) 6(z i-z2)  (8) 

Now  the  sum  and  difference  substitutions 
X  =  (x  !+x  2)  /2 
5  =  X|-x2 

Y  -  (yi+y2)/2 

n  =  yry2 

ks  =  (k  x+k2)/2 

kd  =  k  i~k  2 

and  the  assumption  that  the  frequencies  of  interest  are  much  greater  than 
the  plasma  frequency  so  that 
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01  *  ujj2/  a>!2 

02  =  <^3^/ a)2^ 

enable  Equation  4  to  be  rewritten  as 


where 


—  - - - - [k  y2,  +  1  k  .v2  -  2k  V  .V  ]r 

32  2(k2-k2/4)  d  d  4  d  s  s  s  d 


91/  *♦ 

f-  A(c,n)  - 


k  ,k 


1*  2 


— \ kpA(O) 


r  =  0 


(9) 


V2  = 


7s  *  V 


3S2 

+  _ 

3n2 

32 

♦  il 

3X2 

3Y2 

32 

3  2 

_  + _ 

3^3X  3n3Y 


r  is  now  written  in  terms  of  the  above  sum  and  difference  arguments  as 
rU,n,z,uy)  where  ay  =  ck<j. 

The  unknown  two-frequency  mutual  coherence  function  may  be 
written  as  r  =  r^g  where  r0  is  the  exact  free-space  solution  in  the 
parabolic  approximation 
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which  under  the  sum  and  difference  transformation  used  previously  becomes 

/  i  \2  (  -ik  ( CX+nY)  -  ik  ,[(X2+Y2)/2  +  U2+n2)/8] 

ro  =  _  exp  - - - ° - , - 

lz+zt/  (  (z+zt> 

Substituting  r  =  r jFq  into  Equation  9  and  neglecting  near  zone  terms  of 
the  order  of  ksC  /(z+zj;)  and  smaller,  one  obtains 


ar,  ar, 

+  terms  with  — _  + -  =0  (12) 

3X  3Y 

Equation  12  is  valid  in  the  region  0  £  z  £  t  with  boundary 

condition 

rl(c,n,z=0,«^d)  =  1  (13) 
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Since  the  boundary  condition  is  independent  of  X  and  Y,  and  the  equation 
has  no  terms  other  than  the  derivatives  with  respect  to  X  and  Y,  it  is 
apparent  that 


!Ii  =  ILi  =  o 

3X  3Y 

Now  the  substitutions 

z'  =  2  + 

0  =  C/z'  =  ?/(z+Zt) 
$  =  n/z'  =  n/(z+zt) 
yield  the  following 


3  r  i 

3z  • 


1  k(j  _L  ll_  + 

2  (k^-kj/4)  z'2  L302 


l 


1 

8 


2kjA(z'(0+«)) 
k  ik  2 


The  additional  substitution  =  r2r3  where 

r3  -  exp  j  -  I  A(0)(z'-zt)k^(  L 
(  8  4  p\kl 


(14a) 

(14b) 

(14c) 


(15) 


(16) 


yields  the  simplified  equation  for  r2 


ar2 

3Z  1 


1  k_i  ±_  hi  ♦  •!  lri 

2  z 1  2  [_302  34>2J 


[a( z ■  e,z'  <t>)  -  A(0)  ]r2  =  o 


(17) 


where  kj  has  been  neglected  with  respect  to  ks.  The  effect  of  this 
assumption  is  to  restrict  the  validity  of  the  solution  to  a  small  range  of 
wavelengths  centered  about  ks.  In  the  following  ks,  the  average 
wavenumber,  will  be  replaced  by  k  ,  the  wavenumber  at  the  carrier 
frequency. 


Equation  17  can  only  be  solved  analytically  for  the  case  where 
the  function  A(z'e,z'$)  can  be  represented  as  a  quadratic.  This  is  the 
well  known  quadratic  phase  structure-function  approximation  (Sreenivasiah 
et  al.  (1976);  Sreenivasiah  and  Ishimaru  (1979)).  At  this  point  it  is 
necessary  to  determine  the  form  of  A  for  the  irreqularity  spectrum  of 
interest  here. 

2.2  IONIZATION  IRREGULARITY  DESCRIPTION 

Here  the  ionization  irregularities  are  assumed  to  be  elongated 
along  the  direction  of  the  earth's  magnetic  field.  Figure  2  shows  the 
geometry  of  the  irregularities.  The  electromagnetic  wave  propagates  in 
the  negative  z  direction.  The  magnetic  field  vector  lies  in  the  y-z  plane 
at  an  angle  of  <|*  with  respect  to  the  z  axis.  The  ionization 
irregularities  are  assumed  rotationally  symmetric  with  autocorrelation 
function 


Be(r,s) 


N 


<Ne> 


e 

— 2  exP 


_L !  J 

(qro)2  j 


(18) 
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WAVE 


Figure  2.  A  single  irregularity  elongated  along  the  magnetic  field 
line  in  the  y-z  plane. 

where  s  is  measured  along  the  maqnetic  field  direction  and  r  is 
measured  perpendicular  to  the  magnetic  field  direction.  The  quantity  q 
is  known  as  the  axial  ratio  (Briggs  and  Parkin,  1963)  and  r0  is  the 
correlation  distance  or  irregularity  scale  size. 

The  s,  r  coordinate  system  may  be  related  to  the  x,  y,  z 
system  by  the  equations 

s2  =  (y  sin  ♦  +  z  cos  V»)2  (19) 

r2  =  x2  +  (y  cos  <|/  -  z  sin  *)2  (20) 

Using  the  above  transformat  ion  the  irregularity  correlation  function  may 
be  expressed  in  the  x,  y,  z  system  as 
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This  same  formalism  may  also  be  applied  to  irregularity  power  spectra  that 
are  non-Gaussian.  For  the  case  of  a  power-law  PSD  the  coefficients  Ao 
and  A2  are  different  than  above  but  behave  in  essentially  the  same 
manner  as  a  function  of  the  outer  scale  size  as  long  as  the  three-dimen¬ 
sional  m-situ  electron  density  PSD  falls  off  at  least  as  rapidly  as  K"4. 

2.2.1  Phase  Standard  Deviation 


It  is  useful  to  establish  the  relationship  between  the  coeffi¬ 
cient  Ao  and  the  variance  of  phase  fluctuations  o^2.  For  a  layer  of 
thickness  L  greater  than  the  correlation  length  of  the  irregularities 
the  autocorrelation  function  of  the  phase  is  (Salpeter,  1967) 


<A4>(x1,y1 )  A«|>(x  2  *y2  )> 


Lk2<Ne>2 


4n‘ 


/  B5(Pi-P2,n)  dn 
.00 


Evaluation  of  this  expression  at  pj-p2  =  0  yields 


k£<Ne>2 

4n2 

c 


LAo 


-  knLVk 
.  P  0 


2 

0 


(27) 


(28) 


where  Eguation  6  is  used  for  Aq. 

For  the  anisotropic  correlation  function  given  by  Equation  18,  the 
quantity  Aq  or  A( 0)  is  given  by  Equation  24  and 

a2  =  /¥  (Xre)2q4roLa2  (29) 
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where  the  relationships  nc  =  n/(*2rg)  and  kQ  =  2*/x  have  been  used  in 
Equation  28. 

2.3  ANALYTIC  SOLUTION 

Now  that  expressions  for  Ao,  A2  and  o2  have  been  obtained  it  is 
possible  to  proceed  to  an  analytic  solution  for  r.  Recognizing  that  A(0) 
and  Ao  are  identical,  the  quadratic  structure-function  approximation  may 
be  used  in  Equation  17  to  eliminate  A(0).  With  the  additional 
substitutions 


e  = 


-  k./k2 

2  do 

(30a) 

- 

4  P  0 

(30b) 

(abA2) 1/2  z'  =  a ii' 

(31a) 

(“T — V/4  6  =  a2° 

\a3bA2/ 

(31b) 

(“ F — V/4  ^  =  az* 

\d  bA2/ 

(31c) 

and  the  quadratic  phase  structure-function  approximation  given  by 
Equation  23,  Equation  17  may  be  written  as 


3v  v2j_3p 


r 2  -  ( g2+ Aze2}  v2T2  =  0 


(32) 


With  these  substitutions,  the  boundary  condition  describing  the  mutual 
coherence  function  at  the  transmitter  is 


r2(p,e,v=aizt,u)d)  =  1 


(33) 
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An  analytic  solution  of  the  form 


r2  =  f(v)  exp{-g(  v)u2-h(v)e2}  (34) 

may  be  substituted  into  Equation  32  to  obtain  the  following  three 
equations 

g'  +  +  v2  =  0  (35) 

v2 

h'  +  — —  +  A2v2  =  0  (36) 

v2 

f  +i£f(g+h)_  =  0  (37) 

v2 

Equation  35  represents  the  coefficients  of  terms  from  Equation  32  that  are 
factors  of  u2;  Equation  36  consists  of  the  terms  that  are  factors  of  e2; 
Equation  37  consists  of  terms  that  are  independent  of  g  and  e. 
Additional  simplification  results  from  the  substitutions 


g  =  v2G  (38) 

h  =  v2H  (39) 

which  enable  the  three  equations  to  be  written  as 

G'v2  +  2vG  +  i4 v2G2  =  -  v2  (40) 

H'v2  +  2vH  +  i4 v2H 2  =  -  A2v2  (41) 

f'  +  i  2f ( G+H )  =  0  (42) 
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It  is  not  difficult  to  solve  Equations  40  and  41.  The  results  may  be 
written  as 

G(v)  =  1_  -  1  tanhUv+Cj)  (43) 

4  v  B 

H(v)  =  1-  -  -  tanh(A8v+C2)  (44) 

4v  0 

where 

62  =  -i4  (45) 

and  Cj  and  C2  are  constants  yet  to  be  determined.  With  G  and  H 
known,  f  can  be  determined  by  direct  integration  of  Equation  42  with  the 
result 

f(v)  =  v[cosh(  Bv+C  i)  cosh(A0v+C2)  ]-  */2  C3  (46) 

where  C3  is  another  constant  yet  to  be  determined. 


Equations  38-39  and  43-45  may  be  combined  to  yield  the  result 


for  r, 


r2  =  C3v[cosh(  8v+C  i)  cosh( A3v+C2) ]-1/2 

(  2  ) 
x  exp’-  JL)!  (y2+e2)  +  —  [u2tanh(  Bv+C !)+Ae2tanh( abv+C2)  ]  (47) 

(  4  S  ) 


The  boundary  condition  given  by  Equation  33  is  met  by  requiring 
that  the  coefficients  of  u  and  e  are  zero  at  and  by  requiring 
that  r2  be  unity  at  for  zero  y  and  e.  Mathematically 
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(48) 


- 1  +  tanh(0v.+Ci)  =  0 

4  8 

-iv  A  v2 

-  +  -  tanh(  A0v. +C  2)  =  0 

4  8 

CgV  [cosh( Bvt+C 1)  cosh(Aevt+C2) ]"1/2  =  1 


The  solutions  to  Equations  48  and  49  may  be  found  as 
Ci  =  -  &vt  +  fcn 


w  vty/2 


1  6-1/ v. 


-  8vt  +  Ci 


C9  =  -  ABv*  + 


=  -  A8v^.  +  C2 


AB-l/v. 


(49) 

(50) 


(51) 


(52) 


where  the  new  constants  C\  and  C'2  are  introduced  for  later  con¬ 
venience.  Equations  51  and  52  may  now  be  used  in  50  to  obtain  the 
constant  C3  as 

C3  =  i_  [cosh  C  1  cosh  C’]1/2  (53) 


The  expressions  for  the  constants  Ci,  C2  and  C3  may  now  be  inserted 
into  the  solution  for  r2  with  the  result 
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r2 


.  »  I 


cosh  Cj  cosh  C_2 


I  I/s 


(54) 


( cosh  [e(  v-v^.)+C  i  jcosh  [ab(  v-  )+C 2 ]  j 

1-2  ) 
exp  J  -  —  (u2+e2)  +—  [u2  tanh  [B(  v-v.)+Ci  ]  +  Ae2  tanh [ab( v- v. )+C2  ]  ] i 
(4  3  z  z  ) 


This  result  is  valid  for  propagation  within  the  random  slab  from  z=0  to 
z=L  or  equivalently  from  v=v^.  to  v=vt+v^. 


In  keeping  with  the  approximation  «  k  ,  r2  from  Equation  16 
may  be  written  as 


r3 


I.  A°  ] 

I  *«'  1 


(55) 


At  the  far  edge  of  the  thick  scattering  layer,  z=L  or  z'=zt+L  and  r3 
becomes 


r3  *  exp 


which  may  be  written 


r3  *  exp 


where  o2  is  obtained 
♦ 


A  L  k4  k2  ) 
o  pdf 

8k4  ( 

o  ’ 

as 

•5  “a  I 

H  i 

from  Equation  28. 


(56) 


(57) 
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At  this  point  Equations  54  and  57  may  be  multiplied  to  give  the 
value  of  Tj  at  the  exit  of  the  thick  scattering  layer: 


r^e,  <M'=L+zt,a>d) 


/ L*zt\  | _ cosh  Ci  cosh  Cg _ 1 1/ 2 

\  zt  /  ( cosh ( ea  iL+C i )  cosh( ABC4L+C2) \ 


x  exp 


i(L+zt)a1a2 


(02+<t>2) 


,2  2  2 
+  (L+Zt)  al  a2 
0 


[e2  tanh( 0a  jL+C  1)  +  A$2  tanh( ABa^+Cj)  ] 


! 

I 


(58) 


For  temporary  algebraic  convenience  let  this  be  written  as 
r^e,  <|.,z,=L+zt,a)d)  = 


F  exp {-( A-Bt )  e2  -  (A-ABt ' ) 4>2} 

where 


/ L*zt\  | _ cosh  Ci  cosh  C2 _ ) 

\  zt  /  |  cosh (  Ba xL+C  1 )  cosh(  Aga^+C^)  | 


(59) 


(60) 
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1 


2 


„  i ( L+zt )  aia2 

(61) 

A  4 

,  ,2  2  2 

(L+zt)  aia2 
B  - 

(62) 

0 

t  =  tanh  (3aiL+C 

1) 

(63) 

t'  =  tanh  (AgaiL+C^) 

(64) 

Although  the  above  expressions  are  sufficient,  considerable  simplificat 

occurs  later,  in  the  thin 

phase-screen  approximation,  if  all  the  hyper- 

bolic  functions  of  sums  in 

Equations  60, 

63  and  64  are  expanded.  The 

expansions  require  the  functions  tanh  Ci 

and  tanh  C'2.  These  two 

quantities  are  easily  calculated  as 

tanh  Cl  =  tanh 

i*.  hmH 

1  . 

1 

(65) 

'  Bvt  ' 

»aizt 

f  /A0+l/v  A, 
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tanh  C2  =  tanh 

J 

1 

1 

(66) 

A0aizt 
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The  function  F  can  then  be  written  as 


i  4  ) 

F  =  exp i-  Ba x(L+z. ) 

(  H ) 

x  I _ A _ j1/2 

|  ( 6a  jz.  cosh  8a  jL  +  sinh  Ba  xL)  ( A8a  xz.  cosh  ABa  xL  +  sinh  ABaxL) 

1  (67) 

Similarly 

1  +  8aiZf  tanh  BajL 

t  =  - - -  (68) 

8axZt  +  tanh  8a  iL 


1  +  AeaxZ.  tanh  A8a  xl 

t'  = - - -  (69) 

A8axZt  +  tanh  A6a  xl 


To  complete  the  solution  it  is  necessary  to  solve  for  Tx  in  the 
region  z  >  L  or  z'  _>  L+zt.  Equation  58  serves  as  the  boundary  condition 
at  z'  =  L+z t .  Since  the  region  z'  _>  L+zt  corresponds  to  free-space  with 
no  ionization.  Equation  15  is  appropriate  where  the  last  tern,  which  is 
dependent  on  the  function  A,  may  be  neglected  since  A  is  zero  in  the 
absence  of  ionization  irregularities.  Also,  in  keeping  with  previous 
assumptions,  the  kd2  term  is  ignored  with  respect  to  the  k$2  term  in 
Equation  15. 

The  Fourier  transform  pair 

ri(e,*,z',uid)  =  7  7  el(Ke8+M  ?x(K  ,K  ,z\ud)dK  dK  (70) 

?x(Ke,K  ,z\«>d)  =  ( 1/2 it) 2  7  7  e‘i(K0e+K^)  rx(e,*,z\«  )ded*  (71) 

^  .00  —00 
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may  be  substituted  into  the  suitably  modified  Equation  15  to  obtain  the 
algebraic  equation 


A 

Hi 

3z ' 


+  I 


kd 


2  k; 


1 

2 


(KetK*)ri  ■ 


(72) 


where  k$  has  been  replaced  by  kQ,  the  wavenumber  at  the  carrier  frequency. 

Equation  72  may  be  solved  and  the  boundary  condition  59  applied 
at  z‘  =  L+zt  with  the  result 

riOVVVVV  =  ri(K0,K^L+zt,u»d) 

x  exp  [-i y(K2+K2)  ] 

where 

T .  l  kn  (vL) 

2  ks2  (L*2t)(2tnr) 

The  final  result  may  then  be  obtained  by  taking  the  Fourier  transform  of 

A 

Equation  59  to  obtain  TifKg.K^z'+z^.u^)  and  then  taking  the  inverse 
Fourier  transform  of  Equation  73.  The  appropriate  Fourier  transform  pair 
is  given  by  Equations  70  and  71. 


(73) 

(74) 


The  Fourier  transform  of  Equation  59  is  easily  obtained  as 

fl(K,.Kt.2^*2t>-d)  ’  h 


(A-Bt)1/2  (A-ABt')1/2 


x  exp< 


«S  1  1 

4(  A-Bt )  4(A-ABt')) 


(75) 
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According  to  Equation  73  this  expression  is  multiplied  by  the  appropriate 
factor  to  account  for  free-space  propagation  and  the  inverse  transform 
taken  to  obtain  as  the  final  result 


ri(e,<M'=z t+z  ,» .)  =  - 

[( 1  +  i 4 y ( A -B t ) )  (l+i4y(A-ABt' )) ]1/2 


x  exp 


02(A-Bt) 
1+i 4 y( A-B t ) 


*2(A-ABt')  ) 

l+i4y(A-ABt' )( 


(76) 


Utilizing  the  relationships  given  by  Equations  14b-c  at  the  receiver 
location. 


0  =  c/(zt+zr) 

<J>  =  n/(zt+zr) 

Equation  76  may  be  written  in  final  form  as 


ri(c,n,zt+z  ,«d)  =  - 

[(l+i4y(A-Bt))  (l+i4Y(A-ABt' ))  ]1/2 

j  c2(A-Bt)/(zt+zr)2  n2(A-ABt’ )/(zt+zr)2  | 

XeXPj"  l+i4  y(A-Bt )  ’  1  +  i  4  y  ( A- ABt 1  )  | 

Equation  77  is  the  result  for  after  propagation  through  a  thick  layer. 
The  full  solution  for  the  two-position,  two-frequency  mutual  coherence 
function  is  obtained  by  multiplication  by  r0  as  given  by  Equation  11.  r0 
is  not  affected  by  the  random  layer  and  may  be  ignored  in  the  following. 
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2.4 


PHASE-SCREEN  APPROXIMATION 


Much  simplification  is  possible  if  the  thick  scattering  layer  is 
replaced  by  an  equivalent  thin  phase-screen  of  infinitesimal  thickness  and 
the  same  overall  phase  variance.  To  the  first  order,  in  the  thin  phase- 
screen  approximation 


cosh  Aaa i L  =  1 

(78a) 

Sinh  Aaa  jL  »  Aaa^ 

(78b) 

tanh  Aaa i L  «  AaaLL 

(78c) 

With  the  above  approximat ions  it 

is, easy  to  show  that 

F  |  “d  ) 

pr  2„o2  i 

(79) 

Similarly,  the  coefficient  of  4>2  can  be  written  as 

2 

i(L+zt )ai az 

A-ABt '  * 

4 

A(L+zt)2  aia2  | 

/  2  2  2  V 

(1  +  A  3  a i zt L \ 

3  1 

y  ASa](L+zt)  J 

2  3  2  2 

=  -4  aia2Lzt 

(80) 
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where  e"-  =  -i4  is  used  to  simplify  the  result.  The  quantities  ai  and  a2 
are  qiven  by  Equations  30  and  31  and  may  be  substituted  into  Equation  80 
to  obtain 

A-ABt 1  =  -  A2bA2Lzt2 
A2k4A2Lz2 

=  _  P  t  /Mi 


This  equation  may  be  simplified  by  substitutinq  the  expression  for  o2 

9 

given  by  Equation  28  to  yield 


A-ABt'  = 


A2o2A2z2 


Similarly 


A-Bt  = 


The  final  quantity  of  interest  in  this  phase-screen  approximation  is 


4 y( A-ABt ' )  =  - 


ko  <W  Ao 


L\urh°h 

""o(ztV  flo 


Similarly 


'4 y(A-Bt) 


Wt0^2 

"^tV  Ao 


Now  in  this  phase-screen  approximation  the  result  for  the  two-position, 
two-frequency  mutual  coherence  function  may  be  written  as 


where 


r1(C,n,zt+zr,a)d) 


l  1  1 

/l+i 

fi+i  **  \ 

H  «■) 

i  oj*  /  A2  /  ' 

exp 


l  iiLj  c2/t°  n2/(t°/a2) 

1  vl  \  M  (1+i^)l 


...  <vzr>\ 

0 

’“otztVAo 


■rnV 


Xz..Z4.o.A2 


(86) 


(87) 


(88) 


Equations  86-88  are  valid  for  any  anisotropic  power  spectrum  for  which,  in 
the  stronq  scattering  limit,  the  phase  structure-f unction  can  be  expanded 
in  the  form  of  Equation  23.  For  the  anisotropic  Gaussian  power  spectrum 
of  interest  in  this  work 


Ao 

A2 


(89) 


The  phase  variance  is  given  by  Equation  29  so  that  for  the  Gaussian  PSD 

(90) 


«02  = 


(z.+z  )2r 
t  r  o 


✓ir(  Xre)2qALz2a2 


2»3/2cUt+zr)ro 

r!xVtqsLoN 


(91) 
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2.5 


DECORRELATION  DISTANCE 


Since  is  a  form  of  the  complex  signal  correlation  function  it 
is  possible  to  obtain  the  signal  decorrelation  lengths  in  the  x  and  y 
directions  directly  from  the  exponential  terms  of  Equation  86.  These 
important  signal  parameters  are 


l 

ox 


l 

0 


(92a) 


l 

oy 


(92a) 


40x  and  are  measures  of  the  average  distance  between  fades  at  the 
receiver  location  and  depend  on  the  path  geometry.  The  quantity  A  is 
given  by  Equation  32  and  gives  the  effect  of  the  variation  in  the  angle  of 
propagation  with  respect  to  the  maqnetic  field  direction.  Another  useful 
measure  of  the  signal  decorrelation  properties  is  the  signal  decorrelation 
time  which  is  a  measure  of  the  inverse  fading  rate  or  fading  bandwidth. 

The  signal  decorrelation  time  is  given  by  the  equations 


ox 


oy 


ox 

(93a) 

xeff 

i 

oy 

(93b) 

yeff 


where  ^xeff  and  V  ^  are  the  velocity  components  of  the  1  ine-of-sight  at 
the  receiver  location.  Although  this  velocity  is  a  function  of  transmit¬ 
ter  and  receiver  motion  as  well  as  irregularity  motion,  it  is  sufficient 
to  consider  the  case  where  the  transmitter  and  receiver  are  stationary  and 
the  phase-screen  is  in  motion.  Assume  that  the  scattering  layer  consists 
of  striations  moving  in  unison  at  the  velocity  Vx  i  +  Vy  j  where  i  and  3 
are  unit  vectors  in  the  x  and  y  directions,  respectively.  Consider  only 

the  y  component  of  motion  as  shown  for  thin  screen  case  in  Figure  3.  The 
projection  of  the  striations  at  the  receiver  location  due  to  a  signal 
originating  at  -z^  has  velocity 
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Figure  3.  Thin  phase-screen  propagation  geometry. 


, .  V  (z  +z.) 

V  y  r  t 

yeff  =  -1 - 


(94) 


Using  this  result  for  the  y  component  of  the  velocity  at  the  receiver,  and 
a  similar  result  for  the  x  component  of  velocity,  one  obtains  for  the 
signal  decorrelation  time 


t 

ox 


(95a) 


t 

oy 


(95b) 
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Substituting  the  values  of  Aq ,  A2  and  0^  for  a  Gaussian  PSD  one  obtains  as 
a  final  result 


ox  a  v 
<)>  x 


r  1/2 

0 

it  1  / (qAL )  !/2Xr  V 
e  Ne  x 


(96a) 


oy  ao  v 

<t>  y 


r  !/2 
0 


(%b) 


Tr1/4(qA3L)  1/2Xr  0.  V 
e  Ne  y 

where  both  forms  are  given  for  convenience. 

At  this  point  the  solution  for  the  two-position,  two-frequency 
mutual  coherence  function  is  complete.  The  result  for  the  thick  scatter  - 
inq  layer  given  by  Equation  77  is  simplified  by  takinq  the  thin  phase- 
screen  limit.  In  this  limit  the  result  is  given  by  Equation  86  in  terms 
of  the  quantities  0^,  A,  and  u>' .  These  quantities  are  basic  parameters 
that  describe  the  propagation  environment  and  geometry.  The  quantity  0^ 
is  the  phase  standard  deviation  and  describes  the  strength  of  the  scintil¬ 
lation.  The  decorrelation  distance  of  the  complex  electric  field  is  in 

the  x-direction  and  is  t  /A  in  the  y-direction.  It  will  be  seen  that  the 

0 

parameter  to'  is  related  to  the  coherence  bandwidth  which  is  a  measure  of 
the  maximum  bandwidth  modulated  signal  that  can  propagate  through  the 
layer  without  multipath  distortion. 


In  the  following  the  effects  of  an  aperture  receiving  antenna  on 
the  measured  values  of  signal  decorrelation  distance  and  coherence  band¬ 
width  are  determined. 


ANTENNA  APERTURE  EFFECTS  ON  RECEIVED  SIGNAL 


An  antenna  aperture  acts  to  coherently  collect  the  energy  inci¬ 
dent  upon  the  antenna  and  to  deliver  it  to  the  receiver.  In  the  transmit¬ 
ting  mode  a  directive  antenna  is  designed  to  transmit  energy  only  over  a 
selected  angular  region.  In  the  receiving  mode  this  same  directive 
antenna  accepts  energy  only  from  a  narrow  range  of  angles.  Thus,  relative 
to  an  omnidirectional  point  antenna,  a  directive  antenna  will  experience 
what  is  referred  to  as  angular  scattering  loss  when  the  signal  at  the 
aperture  exhibits  scintillation.  This  angular  scattering  loss  arises  from 
angle-of-arrival  jitter  present  in  the  incident  wavefront  that  may  cause 
energy  to  arrive  at  the  antenna  propagating  at  angles  greater  than  those 
accepted  by  the  receiving  aperture. 

A  different  but  equivalent  way  to  view  the  effect  of  a  receiving 
antenna  aperture  is  as  a  coherent  integrator  of  the  signal  arriving  on  the 
aperture  face.  If  the  antenna  is  pointing  towards  the  source  of  an  undis¬ 
turbed  plane  wave,  then  the  antenna  output  is  maximum.  If  there  are  fluc¬ 
tuations  in  the  signal  phase  across  the  aperture  or  the  incident  wavefront 
is  tilted  relative  to  the  antenna  boresiqht  direction,  the  coherent  inte¬ 
grator  output  is  decreased. 

As  another  aspect  of  the  aperture  averaging  effect,  an  antenna 
can  act  to  increase  the  measured  decorrelation  distance  at  the  antenna 
output.  This  observed  increase  is  caused  by  the  averaging  effect  of  the 
coherent  integration  that  smooths  the  more  rapid  fluctuations. 


Similarly  the  antenna  acts  to  cut  off  signals  that  are  incident 
at  large  angles  from  boresight.  In  severe  scintillation  this  off- 
boresight  energy  propagates  over  longer  paths  than  the  more  direct  energy 
incident  at  boresight.  For  this  reason  a  narrow  beam  antenna  that  does 
not  accept  off-axis  contributions  acts  to  reduce  the  observed  time  delay 
jitter  of  the  output  signal  and  thereby  to  increase  the  observed  coherence 
bandwidth. 

3.1  APERTURE  ANTENNA  FORMULATION 

Consider  the  case  of  a  spherical  wave  that  originates  from  a 
transmitter  located  at  (0,0, -zt)  and  propagates  in  free-space  in  the 
positive  z  direction  where  it  is  incident  on  a  thin  layer  of  irregulari¬ 
ties  located  at  z=0.  After  emerging  from  the  irregularities,  the  wave 
again  propagates  in  free-space  to  a  receiving  aperture  antenna  located  in 
the  plane  z=zr< 

If  the  incident  field  in  the  plane  of  the  antenna  is  U("p,z  ,u>) 
then  the  complex  voltage  envelope  at  the  antenna  output  may  be  expressed 
as  (Price,  Chesnut  and  Burns,  1972) 

v(pQ,zr,w)  =  /  U(T5',zr,u>)  A*(TJ'-TJ0)  d2p'  (97) 

where 

U  =  the  solution  to  the  parabolic  wave  equation 
at  location  T  in  the  receiver  plane  for  a 
monochromatic  signal  frequency  of  w. 

■p  =  the  location  of  the  center  of  the  aperture 
antenna 

it 

A  =  the  complex  antenna  weighting  function 


36 


In  this  section  all  single  and  multiple  integrations  on  infinite  limits 
will  be  written  using  a  single  integral  sign  with  no  limits  shown. 


It  is  recognized  that  Equation  97  is  a  convolution  so  that 


v(K)  =  U(K,zr,u))  A*(K)  (98) 

where  the  Fourier  transforms  are  given  by 

U(K,z  ,u>)  =  ~  /  U(p,z  ,w)  exp(-iK*p)  d2p  (99) 

r  4ir2  r 

A(K)  =  _I_  f  A(p)  exp(-iK.p)  d2P  (100) 

4  TT2 

and  the  vector  K  is 

K  =  k  sine(cos<|>  i  +  sin#  j)  (101) 

for  the  geometry  shown  in  Figure  4. 


The  received  voltage  envelope  at  the  antenna  output  may  then  be 
written  as  the  inverse  Fourier  transform 

v(p0,zr,«)  =  4 it2  /  U(K,zr,u>)  A*(K)  exp(iK*'po)  d2K  (102) 

Since  A(K)  represents  the  antenna  voltage  gain  in  the  direction  K,  it  is 
possible  to  write  the  effect  of  changing  the  antenna  pointing  direction  as 

v(po,zr,u,Ko)  =  4 it2  /  U(K,zr,u>)  A*(K-Ko)  exp(iK*po)  d2K  (103) 
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Figure  4.  Receiving  antenna  aperture  centered  at  origin  of 
coordinate  system. 

where  only  the  antenna  voltage  gain  function  has  been  affected  by  the  new 
pointing  direction.  By  replacing  U  and  A  by  their  Fourier  transforms  and 
performing  the  integrations.  Equation  103  can  be  written  in  the  form 

v(P0,zr,<o,Ko)  =  exp(iKo-po)  /  U(p',zr,w) 

*  A*(p*-po)  exp(-ip'.Ko)  d2p'  (104) 

Equation  104  gives  the  voltage  output  from  an  aperture  antenna  centered  at 

p  and  pointed  in  the  K  direction.  This  result  is  applicable  to  the  case 
o  o 

of  a  transmitted  monochromatic  sinusoidal  waveform. 

Now  if  the  transmitted  waveform  is  a  pulse  waveform  modulated  on 
a  carrier,  the  transmitted  signal  can  be  expressed  as 
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(105) 


s ( x)  =  Re  {m(  t)  exp(  i  oiq  x)  } 

where  o>0  is  the  carrier  radian  frequency  and  m(  t)  is  the  modulation  wave¬ 
form.  After  passage  through  a  layer  of  irregularities  the  received  signal 
may  be  written  (Knepp,  1982)  in  the  absence  of  an  antenna,  as 

r(-p,zr,T)  =  Reje(p,zr,x-Td)exp[i(wot+eo}J  (106) 

where 

1  00 

e(p,z  ,t)  = —  /  M(w)U(  p,z  ,w<-u)  )exp(iwT)dw  (107) 

r  2n_. 

The  quantity  e  is  called  the  complex  envelope  of  the  received  waveform. 
M(w)  is  the  Fourier  transform  of  the  transmitted  modulation  waveform  and 
is  given  by 


M(u>)  =  /  m(  T)exp(-i«)T)d  t  (108) 

—  00 

The  quantities  90  and  are  the  mean  phase  shift  and  time  delay 
incurred  after  propagation  through  the  ionization  and  are  given  as 

e i  «  -  JSL  /r(l-u>W)i/2  dz.  (109) 

c  -zt 


U!2/  ID2) 

p  o' 


l/2 


dz' 


(110) 
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In  order  to  obtain  the  effect  of  the  antenna  on  the  received 
time-domain  waveform,  the  received  complex  envelope  L)  in  Equation  107  is 
replaced  by  its  equivalent  value  v  given  by  Equation  104.  Thus  the 
received  time-domain  signal  at  the  antenna  output  may  be  written  as 

e(p  ,z  ,t,K  )  =  I -  /  M( ui)U('pl  ,z  , o>+o)  ) 

U  i  (J  o  i  U 

2ir 

x  A  ("p1 -7^) exp (-i’p'  •Ko+iwT)d2p'du)  (111) 

where  the  phase  term  exp(iK  *p  )  has  been  omitted  since  it  is  constant. 

oo 

In  the  case  that  the  receiving  aperture  is  a  point  antenna.  A*  is  a  delta 
function  and  the  received  complex  envelope  given  by  Equation  111  is 
identical  to  the  result  given  by  Equation  107. 


The  correlation  function  of  the  received  power  is  given  from 
Equation  111  as 

<e(p1,T)e*(pz,T)>  =  L  /  M(to'  )M*(u)")<U(’pl,aj'+(i)  )U*(p,,,«,,+m J> 

4  IT2  0  0 


x  A  ( p1 -Pi)A( p" - p2) 


x  exp  [-i(  p'-p")  »K  +  i ( u>' -a>" )  T]du>'dw"d2p,d2p" 


(112) 


where  the  zr  dependence  has  been  omitted  and  it  is  assumed  that 
interest  is  focused  on  only  the  K0  pointing  direction.  Recognize  that 
the  term  in  the  angle  brackets  is  the  two-position,  two-frequency  mutual 
coherence  function,  and  is  only  a  function  of  the  difference  in  position 
and  frequency,  rCp'-p’\  <i>'-ui").  With  this  in  mind  introduce  the  sum  and 
difference  transformations 
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(113) 


0),  =  Oi'  -  <ii" 

(114) 

UJ  =  —  (  UJl  +  {j)"  ) 

S  2 

to  obtain  for  the  correlation  function  of  received  power 

<e(p!,  t)e*(  P2,  t)>  =  /  r(pd,wd)M(a)s+u)d/2)M  (u>s-u>d/2) 

x  A  (ps+pd/2-p1)A(ps-pd/2-p2) 
x  exp( -i"pd  *K0+i  u>d  T)dwddu)5d2pdd2ps 

(115) 

In  order  to  simplify  this  equation  it  is  convenient  to  consider 
the  case  that  the  transmitted  power  is  a  delta  function  of  delay.  Since 
Vl(io)  is  the  Fourier  transform  of  the  transmitted  complex  envelope,  the 
time-domain  representation  is  given  by 

m(  x)  =  —  /  M(u))exp(  i  wx)d(u  (116) 

2  w 


Thus  the  input  power  can  be  written  as 


--  n(t)ri*(t) 


=  -i—  /  M(  w1  )M(  u>" )exp  [i  ( ui1  -<i>"  )t  Jdw'dw" 


(117) 
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Now  under  the  sum  and  difference  transformations  above  the  input  power  can 
be  written  as 


I(t) 


4tt< 


/  M(u>s+u>d/2)M 


( V  “d^  exp  ^ 1  wd  * ) d  “dd  “s 


(118) 


Now  if  the  input  power  is  a  delta  function 

I(t)  =  <$(t)  =  L  /  exp(  i  ojdt)do»d  (119) 

2  IT 

Comparison  of  Equations  118  and  119  yields  for  an  input  power  delta 
function, 

—  /  M(o)  +oj,/2)M  (  id  —  oj . / 2 ) d a>  =  1  (120) 

2s  a  so  s 

IT 


Equation  120  above  may  be  used  to  obtain  the  received  power 
correlation  function  for  an  input  power  delta  function 

G(‘p1-'p2»t)  =  —  / 

2  it  d  a 

x  a  ( Ps+Pd/2-^i)A('ps-'pd/2-'p2) 

x  exp(  -i“pd  *^0+1u,dT)dudd2Pdd2Ps  (121) 

where  the  notation  G  is  used  for  this  important  function.  At  this  point 
introduce  the  Fourier  transform  of  the  antenna  aperture  distribution  as 
given  by  the  Fourier  transform  of  Equation  100 

A(p)  =  /  A(K)  exp  (iK*p)d2K  (122) 

Substitution  of  Equation  122  into  Equation  121  gives 
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G(?i-P2,t)  =  L  /  r(pd,a.d)A*(K1)A(K“) 


2* 


exp[-iP  .(K'-K")-i  fl.(K'+K") 
s  2 


+  iK'  •Pi-iK"  •  P2 - "i pd ”Ko+1  tJ°dT ] d^fjd 2 Pc)daPsd2K  'd2K" 

(123) 

The  ps  integration  yields  the  delta  function  4it26(k  1  -K" )  and  enables 
one  to  then  easily  perform  the  K"  integration.  The  result  is  then 

G(?i-U2,t)  =  2*  /  r(T5d,a>d)|A(K)|2 


x  exp  [iK*(Pi-P2)-i(K+Ko)*Pd+iwdT]dwdd2Pdd2K  (124) 

To  additionally  simplify  this  expression  introduce  the  generalized  power 
spectrum  (Wittwer,  1980;  Knepp  1982)  which  is  the  Fourier  transform  of  the 
mutual  coherence  function, 

S(K,t)  =  _L  /  r(TTd,«d)  exp  (-iK-Pd+iu)dT)da)dd2Pd  (125) 

8n 

Utilization  of  the  generalized  power  spectrum  in  Equation  124  yields  for 
the  power  impulse  response  function  at  the  aperture  output 

G(pi-p2,t)  =  (2»)4  /  S(K+K',t)|A(K)|2  (126) 

x  exp  [  iK*  ("Pi  -7J2 )  ]d2K 


A  simple  change  of  variables  yields  the  final  result 
G(?i-TS2,t)  =  (2»)4  /  S(K,t)|A(K-K0)|2 

x  exp  [i(K-K0)*(Pi-p2) ]d2K  (127) 
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Taking  Pi=P2  gives  the  mean  power  at  the  antenna  output  due  to  a  trans¬ 
mitted  power  delta  function 

G(0,t)  =  (2w)4  /  S(K,T)|A(K-KQ)|2d2K  (128) 

It  is  easy  to  verify  that  the  above  equation  has  the  correct  normaliza¬ 
tion.  For  a  point  receiving  antenna  ACp)  =  6fp)  and  therefore 

A(K)  «  -L  (129) 

4ir2 

and  the  mean  received  power  impulse  response  function  becomes 

G(0, t)  =  /  S(K,t)d2K  (130) 

which  is  identical  to  previous  results  (Knepp,  1982)  that  were  derived  for 
the  case  of  no  antenna  or  equivalently,  a  point  antenna. 

Equation  128  can  be  written  in  another  very  useful  form  with  the 
substi tution 


r(K,wd)  =  ~  /  r(pd,<i)d)exp(-iK*pd)d2pd  (131) 

4  IT 

A  __ 

where  r(K,oy)  is  the  Fourier  transform  of  the  mutual  coherence  func- 

A 

tion.  If  r  is  used  in  Equation  124  one  obtains 

G(0,t)  «  (2 ir) 3  /  r(K+K0,ood)jA(K)' 2exp(iu)dT)du)dd2K  (132) 
or,  changing  variables, 

G(0,  t)  »  (2w)3  /  r(K,Wd)|A(K-K0)|2exp(iUdT)da,dd2K  (133) 

It  will  be  seen  that  Equation  133  is  preferable  for  the  evaluation  of  the 
effects  of  antennas  on  mean  time  delay  and  time  delay  jitter. 
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It  is  easy  to  compute  f  directly  from  Equation  131  using  the 
thin  screen  approximation  for  the  two-position,  two-frequency  mutual 
coherence  function  given  by  Equation  86.  The  result  is 


r(K,u»d)  =  —  exp 


4itA 


l.fMl 

I  2“o  ) 


exp 


(  k2*7?  /  U)h\  /  0),  \) 

-  JL2.  i+i  _!  .  jjl  (i+i 

I  4  \  w1  /  4a2  \  w'/A2/| 


(134) 


The  generalized  power  spectrun  may  now  be  computed  by  using  the  expression 
just  given  for  f  and  by  applying  the  relationships  of  Equation  125  and 
131  to  obtain 


where 


a  = 


a  u‘ 


(135) 


(136) 


For  a  Gaussian  irregularity  power  spectrum,  the  parameter  a  may  be 
calculated  using  the  expression  for  u>'  given  by  Equation  88  as 


a 


xVt% 

"<Vzr)r 


2 

0 


(137) 


It  is  apparent  that  a  is  a  measure  of  the  severity  of  the  scattering 
since  it  is  directly  proportional  to  the  phase  standard  deviation. 


3.2 


APERTURE  DISTRIBUTION  FUNCTION 


Let  the  antenna  be  modeled  as  a  Gaussian  beam  with  gain  function 
Ga(e)  =  Go  exp  (-e2/e2)  (138) 

where  the  beamwidth  90  is  related  to  the  effective  aperture  diameter 
D  by 

2 

9o  =  (^)/(4  *n2)  (139) 

is  the  wavelength  at  the  carrier  frequency.  With  this  choice  of 

Ga{V2D)=Go/2  (140) 

so  that  this  Gaussian  aperture  has  approximately  the  same  3  dB  beamwidth 
as  a  uniformly  illuminated  circular  antenna. 

The  antenna  gain  pattern  is  the  square  of  the  transform  of  the 
aperture  response  function  defined  earlier.  Thus 

A(K)  =  (Ga( 0)  |  =  G^2  exp  (-e2/2  02)  (141) 

Equation  122  may  be  used  to  find  the  aperture  weighting  function  A(T>). 
Then  since 


where  ° 


K*'?  =  k  sin6(x  cost  +  y  sint) 


(142) 


dK  =  k2  sinQ  cosQ  d®dt 

the  substitutions 


(143) 


X  =  P  cost' 
y  =  p  sint1 
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yield,  after  performing  the  azimuthal  integration 

TT 

A("p)  =  2*k2G*/2  /  Jo(kPsin9)exp(-02/202)  sin0  cos0  d0 
o 

For  small  0O  most  of  the  contribution  to  this  integral  comes  from 
small  values  of  9.  Thus  it  is  possible  to  replace  the  sine  of  the  angle 
with  the  angle  itself  and  to  extend  the  upper  limit  to  infinity 
(Ratcliffe,  1956).  The  resulting  expression  is  available  in  standard 
integral  tables  (Gradshteyn  and  Ryzhik,  1965)  with  the  result 

A(p)  =  2*k262Gj/2  exp  (-kzp202/2)  (145) 

Since,  in  the  small  angle  approximation, 

K2  =  k2  sin2  0  «  k202  (146) 

Equation  141  may  be  written  as 

A(K)  =  G^/2  exp  (-Kz/2kz0z)  (147) 

It  is  easy  to  verify  that  Equations  145  and  147  are  Fourier  transform 
pairs . 

3.3  ANGULAR  SCATTERING  LOSS 

The  scattering  loss  and  the  effect  of  the  antenna  on  measure¬ 
ments  of  decorrelation  distance  may  both  be  calculated  simultaneously  by 
evaluation  of  Equation  127.  Here  only  the  limit  of  strong  scattering 
(large  o<$,)  and  large  coherence  bandwidth  (large  u>coh  or  u>’)  is  consider¬ 
ed.  In  that  case  the  parameter  a  given  by  Equation  136  is  small  and  the 
generalized  power  spectrum  reduces  to  the  simplified  form 
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1  im  S(K,t) 

a+O 


*2otu)'  (  K2*2  K 2l2  ) 

_  exD  L  x  0  -  y  0  ! 

25/2*3/24  P(  4  4A2  ) 


x  exp 


\  -  -  (aw' t)2 
(  2 


I 

I 


(148) 


The  factor  au>'  is  equal  to  the  ratio  0^/04,  and  does  not  reduce  to 
zero  in  the  limit  of  small  a. 


The  correlation  function  of  the  received  power  can  easily  be 
evaluated  for  the  case  of  a  Gaussian  antenna  beam.  Using  Equations  147 
and  148  in  Equation  127  one  obtains 


( 2 * )  4G  <W  exP  T> 

G(pi-p2»t)  =  - - - ? - 

(2ir)  l/2  [(1+0. 28D2/ Jl2)  (1+0.28a2D2/ l2)  l1/2 

(  -42/2o  n2/(i2/A2)  j 

X  exp  - - — - - — - 

l+0.280z/£2  1+0.28a2D2/ A2) 

'  O  0 


(149) 


where  it  is  assumed  that  the  antenna  boresight  is  in  the  z-direction  and 
therefore  Ko=0.  In  Equation  149 

"Pi  -  7>2  =  (4»n) 

and  the  antenna  beamwidth  e0  is  written  in  terms  of  x/D  as  given  by 
Equation  139.  For  an  omnidirectional  receiving  antenna  whose  aperture 
diameter  approaches  zero,  the  gain  G0  goes  to  l/(2n)4  and 

lim  6(0, t)  -  - exp  j  -  I  (du.'T)2  |  (150) 

(2w)l/2  I  2  I 


which  is  identical  to  earlier  results  for  the  power  impulse  response 
function  neglecting  antenna  aperture  effects  (Knepp,  1982). 
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The  ratio  of  the  correlation  function  of  the  received  signal 
envelope  (Equation  149)  to  the  received  power  with  no  antenna  (that  is, 
with  a  wide-beamwidth  point  antenna)  gives  the  aperture  antenna  effect  as 


G(p!-P2»t)  _  _ ^ _ 

lim  G(0,  t)  [(1+0.28D2/ i.^)(l+0.28A2D2/t^)  j1/2 

D-0 

(  ;2/*2  n2/(*2/A2)  | 

X  6XP  j  1+0. 28D2/ l2  1+0. 28A2D2/ i2 1 


(151) 


The  angular  scattering  loss  is  given  by  the  term  in  the  square  bracket 


loss  =  [(1+0.28D2/A2)(1+0.28a2D2/*2) ]1/2  (152) 

Equation  152  gives  the  angular  scattering  loss  as  a  function  of  antenna 
diameter  and  various  properties  of  the  scattering  medium  and  scattering 
geometry.  In  order  to  separate  geometric  effects  from  antenna  effects 
note  that  the  decorrel ation  distance  (with  no  antenna)  is  given  by 
Equation  90  repeated  here 

(z.+z  )2r 

l2  =  t  r;  o 

°  ✓*  Urg)2qALz2o2 

where  the  factor  A  gives  the  effect  of  varying  the  inclination  angle. 
For  propagation  parallel  to  the  magnetic  field  line  «|f0  and  A-l. 
Therefore  define  the  new  parameter  as  the  decorrel ation  distance 

for  parallel  propagation  as 


(zt+zr}  o 
(Arg)2qLz2o2 


(153) 
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As  can  be  seen  ip  is  not  a  function  of  the  propagation  geometry.  Now 
the  angular  scattering  loss  can  be  written  in  terms  of  £p  as 

loss  =  [(l+0.28AD2/**)(l+0.28A3D2/iJ)]1/2  (154) 

which  explicitly  gives  the  effect  of  variations  in  the  antenna  aperture 
size  and  of  variations  in  the  scattering  geometry. 


As  seen  from  Equation  151  and  the  following,  the  effect  of  the 
antenna  aperture  is  easily  interpreted  as  a  function  of  the  ratio  of  the 
antenna  diameter  to  the  decorrelation  distance  of  the  field  incident  on 
the  aperture.  It  is  easily  shown  that  these  results  may  eguivalently  be 
interpreted  in  terms  of  the  ratio  of  the  rms  scattering  angle  to  the 
antenna  heamwidth.  From  the  Appendix  and  Equation  139  relatinq  the 
antenna  beamwidth  and  diameter  one  obtains 


D 

7T 

°0x 

l 

0 

(21n2) 

1/2 

0 

0 

AD 

7T 

% 

l 

0 

(21n2) 

1/2 

0 

0 

where 

°0x 

and 

OQy  are 

the 

rms  angle 

-of-arrival 

along 

the 

X- 

and  y-axes. 

respective 

iy. 

Thus 

since 

0.28  D2  =  2 

■5, 

l2  02 

o  o 


(155) 

(156) 

fluctuations  measured 

(157) 
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(158) 


0.28  A2D2  _  2  CT0y 


it  is  apparent  that  antenna  aperture  effects  become  important  only  when 
the  rms  scattering  angle  approaches  or  exceeds  the  antenna  beamwidth. 

This  statement  is  true  for  all  the  effects  considered  here;  the  severity 
of  the  effect  depends  upon  the  particular  measurement. 

Figure  5  shows  the  scattering  loss  in  decibels  as  a  function  of 
the  field  inclination  angle  for  a  value  of  the  axial  ratio  q  of  15.  In 
all  succeeding  figures  q  is  taken  as  15  (Wittwer,  1979).  The  ratio  of 
the  antenna  diameter  to  the  decorrel ation  distance  for  parallel 
propagation  D/tp  is  shown  parametr ical ly  for  the  values  1,  3,  10,  30 
and  100.  Figure  6  is  another  plot  of  Equation  154  and  shows  the  angular 
scattering  loss  as  a  function  of  the  normalized  antenna  diameter  D/ ip 
for  values  of  the  inclination  angle  of  0°,  15”,  30°,  45°  and  90°.  It  is 
seen  in  both  figures  that  only  when  the  antenna  diameter  is  large  with 
respect  to  the  decorrelation  distance  is  the  angular  scattering  loss 
significant.  Small  inclination  angles  cause  increased  values  of  the  phase 
standard  deviation  and  thus  give  increased  angular  scattering  and 

increased  scattering  loss  as  shown.  An  increase  in  the  inclination  angle 
causes  a  decrease  in  o$,  an  effective  increase  in  the  true 
decorrelation  distance  and  therefore  causes  a  decrease  in  the  effect  of 
angular  scattering  loss. 
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3.4 


SIGNAL  DECORRELATION  DISTANCE 


The  effect  of  the  aperture  on  the  measured  correlation  function 
is  given  by  the  two  terms  in  the  exponential  of  Equation  151  as 

*x  =  yi+0.2802/*2)  V2  (159) 

l  =  ( t  /A)(l+0.28A2D2/i2)  V2  (160) 

y  o  o 

where  ix  is  the  antenna  measured  decorrel ation  distance  in  the 
x-direction  and  iy  is  the  same  quantity  measured  in  the  y-direction. 
These  quantities  differ  because  of  the  geometry  of  propagation  with 
respect  to  the  magnetic  field  direction.  For  propagation  parallel  to  the 
magnetic  field  line  A=1  and  £x=4y.  Similarly  for  isotropic  irregu¬ 
larities  q=l  and  A=1  and  0  is  not  a  function  of  the  field  inclina¬ 
tion  angle  so  that  zx=ty  again.  With  a  point  receiver,  D=0,  and  the 
measured  decorrel at  ion  distances  are  Iq  and  £q/a  as  given  by 
Equations  92a  and  92b  as  well  as  Equations  159  and  160. 


The  decorrelation  distances  in  the  x-  and  y-directions  may  be 
normalized  by  their  values  measured  by  an  omnidirectional  (point)  antenna 
and  plotted  as 


—  =  (1+0.28AD2/*2) l/2 
lo  P 

-J?—  =  ( 1+0. 28A3D2/ L2)  */2 

(VA)  P 


(161) 

(162) 


Figure  7  shows  the  relative  decorrelation  distance  in  the  x-direction  as  a 
function  of  inclination  angle  for  values  of  D/jtp  of  1,  3,  10,  30  and 
100.  The  axial  ratio  q  is  taken  as  15.  Figure  8  shows  the  relative  de- 
correlation  distance  in  the  x-direction  as  a  function  of  the  normalized 


RELATIVE  DECORRELATION  DISTANCE  (X) 


antenna  diameter  for  values  of  the  inclination  angle  of  0*,  15*,  30*,  45* 
and  90*.  Fiqure  8  shows  that  antennas  with  diameters  less  than  the 
decorrelation  distance  do  not  affect  measurements  of  anqular  scattering. 
From  Figures  7  and  8  it  is  evident  that  very  large  antennas  can  qreatly 
increase  the  measured  decorrelat ion  distance  (relative  to  a  point  antenna) 
for  sufficiently  strong  anqular  scattering  at  small  inclination  anqles. 

Figures  9  and  10  show  the  effect  of  an  aperture  antenna  on 
measured  values  of  the  decorrelation  distance  in  the  y-direction.  The 
format  of  these  two  figures  is  identical  to  the  format  of  Figures  7  and 
8.  In  the  y-direction  the  effect  of  the  antenna  aperture  is  much  smaller 
than  in  the  x-direction  for  values  of  inclination  angle  greater  than  about 
10*. 


3.5  ANTENNA  APERTURE  EFFECT  ON  <t>  and  ox 


In  the  same  manner  than  an  aperture  antenna  with  a  finite 
beamwidth  neglects  or  averages  out  enerqy  incident  at  of f -bores iqht  anqles 
to  reduce  the  received  power,  it  also  acts  to  reduce  the  measured  time 
delay  and  time  delay  jitter.  This  reduction  occurs  because  the  energy 
arriving  from  directions  away  from  boresiqht  typically  travels  over  lonqer 
paths  than  the  more  direct  signal.  These  lonqer  paths  require  a  lonqer 
propagation  time  and  hence  contribute  to  increased  <x>  and  oT  values. 

If  this  enerqy  is  neglected  by  an  aperture  antenna,  then  the  siqnal  at  the 
output  will  be  character ized  by  smaller  <t>  and  oT  than  that  measured  by 
an  omnidirectional  (point)  antenna.  ^ 


Now  define  the  mean  time  delay,  <t>,  as 

<t>  =  /  G(0»T)Tdx 
I  G(0,T)dx 


(163) 
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Figure  9.  Effect  of  aperture  and  magnetic  field  direction  on  measurements 
of  the  decorrelation  distance  in  the  y-direction. 


Figure  10.  Effect  of  aperture  and  magnetic  field  direction  on  measurements 
of  the  decorrelation  distance  in  the  y-direction. 


and  define  the  time  delay  jitter  oT,  by  the  second  moment 

02  ,  /  G(0,T)T2dr  _  <t>2  (164) 

j  G(0,x)dT 

where  G(0,t)  is  the  response  measured  at  the  output  of  the  receiving 
antenna  due  to  an  impulse  of  power  originating  from  the  transmitter.  Now 
the  alternative  expression  for  G(0,t)  given  by  Equation  133  is  the  choice 
for  the  evaluation  of  the  mean  time  delay  and  the  time  delay  jitter.  In 
this  case  the  integrals  with  respect  to  r  give  Dirac  delta  functions  and 
derivatives  of  Dirac  delta  functions.  That  is 


J  r(K,a>d) 

exp  ( 

iu>dT)dx  =  2irr(K,0) 

(165) 

/  r(K,cod) 

exp  1 

[  i  u>d  t)  rdt  =  i2n 

ar(K,u) 

3u.d 

O 

It 

3° 

(166) 

/  r(K,a>d) 

exp  1 

[  i  ud  t)  T2d  T  =  -  2lT 

32r(K,u^) 

3a» .  _n 

d  u) .  =0 

(167) 

Utilizing  Equations  133  and  165-167  in  the  expressions  for  the  mean  time 
delay  and  the  time  delay  jitter  leaves  only  integrals  over  K.  For  the 
Gaussian  antenna  beam  used  previously  these  integrals  are  easily  performed 
with  the  result  that,  in  terms  of  the  antenna  diameter 


<t>  = 


_ ♦ _ *? _ I 

2u’  |1+0.28D2/*2  1+0.28A2D2/*2 j 


(168) 


and 


1  +  M  1  +  )  1 

a2u»'  2  2  ((1+0.28D2/ i2)  2  (1+0.28A2D2/*2)  2  j  <*>'  2 


(169) 
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Since  awl=o^/u^  the  first  term  in  Equation  169  is  proportional  to 
the  ratio  of  the  phase  standard  deviation  and  the  carrier  radian  fre¬ 
quency.  This  term  is  always  small  for  GHz  frequencies  and  may  be 
neglected.  The  result  for  the  time  delay  jitter  is  then 


1 


1 


a/2 


T  — 


/2  o'  I (1+0.28D 2/£2)2  (1+0.28A2D2/£2)2I 


(170) 


For  point  antennas  with  no  directivity  the  mean  time  delay  is  given  from 
Equation  170  as 


<x(point  antenna)>  =  —5^  (1+A2) 

2w' 


(171) 


Thus  the  effect  of  the  aperture  antenna  on  the  measured  mean  time  delay  is 
given  by  the  ratio 


<x(with  antenna)> 
<x(point  antenna)> 


(1+0.28A2D2/*2)  +  A2( 1+0. 28D2/ £2) 

' _ o' _ _ o 

(1+A2)(1+0.28D2/ A2) ( 1+0. 28a2D2/ £2) 


(172) 


For  a  point  antenna  the  time  delay  jitter  is  given  by  Equation  170  as 


oT(point  antenna)  =  ^ —  (1+A4)1/2  (173) 

/2  u)' 


The  effect  of  the  receiving  aperture  is  best  measured  by  the  ratio 


°x(with  antenna) >  _ 
o^(point  antenna)> 


f ( 1+0 . 28 A2D2/ £2) 2  +  A4( 1+0.28D2/*2) 2}1/2 
(1+A4) 172( 1+0.28D2/ t2)(l+0.28A2D2/t2) 


(174) 
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3.5.1 


Coherence  Bandwidth 


The  coherence  bandwidth  is  the  inverse  of  the  time  delay  jitter 
of  the  power  impulse  response  function.  From  Equation  170 


n  u>* ( 1+0 . 28D2/ i2 ) ( 1+0 . 28A2D 2 /l2 ) 
_ o_ _ o 

{( 1+0 . 28A2D2/ t2) 2  +  A4(1+0.28D2/*2)2}1/2 


(175) 


With  a  point  antenna,  the  measured  coherence  bandwidth  becomes 


«cot,(poi"t 


antenna)  = 


/2  w' 

(1+aV/2 


(176) 


For  isotropic  irregularit ies,  q=l  and  therefore  A=1  so  that  “coh^' 
as  is  given  by  Knepp  (1982).  For  infinitely  elongated  irregularities  q 
goes  to  infinity,  A  tends  to  zero  and  a)co^='^  which  also  agrees 
with  the  previous  results  that  were  calculated  without  consideration  of 
aperture  averaging. 


The  effect  of  an  aperture  antenna  on  the  measured  coherence 
bandwidth  is  given  by  the  ratio  of  the  antenna  measured  coherence  band¬ 
width  to  the  coherence  bandwidth  with  no  antenna.  This  ratio  is  the 
reciprocal  of  Equation  174.  Figures  11  and  12  show  the  effect  of  inclina 
tion  angle  and  of  aperture  diameter  on  the  measured  coherence  bandwidth. 
In  both  figures  the  ordinate  is  the  relative  coherence  bandwidth  given  by 
the  reciprocal  of  Equation  174.  In  Figure  11  the  relative  coherence  band 
width  is  shown  as  a  function  of  the  inclination  angle  between  the  direc¬ 
tion  of  propagation  and  the  magnetic  field.  Values  of  the  ratio  of  the 
antenna  diameter  to  the  decorrelation  distance  for  parallel  propagation 
are  shown  parametrically.  In  Figure  12  the  abscissa  and  the  parametric 
quantity  are  exchanged.  The  two  figures  show  that  the  aperture  antenna 
has  a  strong  effect  in  increasing  the  measured  coherence  bandwidth  over 
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the  value  given  by  an  omnidirectional  antenna  especially  for  large  aper¬ 
tures  and  propagation  angles  close  to  the  magnetic  field  direction.  This 
increased  measurement  of  coherence  bandwidth  is  caused  by  the  action  of 
the  aperture  to  exclude  or  average  out  signals  that  are  incident  at  large 
angles  and  hence  have  experienced  greater  delay  than  the  directly  incident 
signals.  Equivalently,  as  discussed  in  conjunction  with  Equations  155- 
158,  it  is  apparent  that  the  aperture  antenna  significantly  affects  mea¬ 
surement  of  the  coherence  bandwidth  whenever  the  rms  angle-of-arr ival 
fluctuation  is  greater  than  the  antenna  beamwidth. 

Note  that  Figures  7-12  show  the  measured  value  of  some  signal 
parameter  (measured  with  an  aperture  antenna)  normalized  by  the  value  of 
that  parameter  measured  by  an  omnidirectional  antenna.  That  is,  these 
figures  present  the  effect  of  the  aperture.  Actual  values  of  these  para¬ 
meters  are  easily  calculated  for  any  ionospheric  condition  and  propagation 
geometry  using  the  equations  provided. 
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SECTION  4 
SUMMARY 


In  this  report  results  are  presented  that  show  the  effect  of 
non-isotropic  irregularities  on  the  two-position,  two-frequency  mutual 
coherence  function.  In  addition  the  aperture  averaging  effect  of  a 
Gaussian  antenna  is  determined  for  strong  scattering  conditions.  It  is 
shown  that  antennas  that  are  larger  in  diameter  than  the  decorrelation 
distance  can  experience  significant  angular  scattering  loss  and  cause 
increased  measurements  of  decorrelation  distance  and  coherence  bandwidth 
relative  to  values  measured  with  an  omnidirectional  antenna.  Increases  in 
measured  decorrelation  distance  are  caused  by  the  averaging  effect  of  the 
aperture  that  smooths  the  smaller  scale  signal  fluctuations.  Increased 
measurements  of  coherence  bandwidth  are  caused  by  the  action  of  the 
antenna  to  cut  off  signal  contributions  originating  from  off-boresight  and 
which  generally  have  experienced  more  delay  than  the  direct  signal  inci¬ 
dent  at  bores ight. 

Generally,  the  effects  of  aperture  averaging  begin  to  be  signif¬ 
icant  for  aperture  sizes  that  are  ten  times  the  decorrelation  distance. 

For  still  larger  aperture  dimensions,  the  aperture  has  a  significant 
effect  on  the  parameters  that  describe  the  received  signal  at  the  antenna 
output.  However,  this  effect  may  be  calculated  using  the  formulas  and 
curves  provided  herein. 
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APPENDIX 

EVALUATION  OF  RMS  ANGLE-OF-ARR I VAL  FLUCTUATION 


Consider  a  plane  wave  traveling  in  the  negative  z-direction  and 
incident  in  the  x-y  plane  at  an  angle  0  from  the  z-axis.  The  electric 
field  is  given  by 

E(x,y)  =  Eq  exp  fik(sin0  x+cos0  z)}  (A-l) 

where  the  exp(iwt)  time  dependerice  is  suppressed.  The  anqle-of-arrival 
measured  along  the  x-axis  is  computed  as 


e  =  1 

9E(x,z) 

*  ikE 

3x 

0 

This  computation  gives  sin  0  which  in 

approximately  9. 

In  the  case  of  interest  here 

solution  to  the  parabolic  wave  equation 

measured  in  both  x- 

and  y-directions  as 

e  =  * 

3U{?,Z,(*») 

x  ikUQ 

3x 

0  =  1 

3U(?,z,w) 

y  ikUQ 

3y 

(A-2) 


(A-3) 


(A-4) 


fi5 


where  the  meausrement  depends  on  the  frequency.  At  a  single  frequency  the 
rms  values  of  9X  and  0 y  are  related  to  the  second  derivatives  of  the 
two-position  mutual  coherence  function  (Papoulis,  1965,  p.  317) 


1  32r(c,n,u)d=0) 

k1  3C1 


l  a2r(c,n,u.d=o) 

k2  an2 


(A-5) 


(A-6) 


From  Equation  86  in  the  text  one  obtains  for  the  rms  angle-of-arrival 
fluctuations 


o 


0x 


/2 

k2o 


(A-7) 


/2 

k(aQ/A) 


(A-8) 
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ATTN:  CCC-EMEO-PED,  G.  Lane 

ATTN 

FC-1 

US 

Army  Communications  Cmd 

ATTN:  CC-OPS-WR,  H.  Wilson 

Field  Command 

Defense  Nuclear  Agency 

ATTN:  CC-OPS-W 

ATTN 

FCTT,  W.  Summa 

US 

Army  Communications  R&D  Cmd 

ATTN 

FCTXE 

ATTN:  DRDCO-COM-RY,  W.  Kesselman 

ATTN 

FCTT,  G.  Ganong 

ATTN 

FCPR 

US 

Army  Foreign  Science  &  Tech  Ctr 

ATTN:  DRXST-SD 

Interservice  Nuclear  Weapons  School 

ATTN 

TTV 

US 

Army  Materiel  Dev  *  Readiness  Cmd 

ATTN:  DRCLDC,  J.  Bender 

Joint  Chiefs  of  Staff 

ATTN 

C3S 

US 

Army  Nuclear  &  Chem  Agency 

ATTN 

C3S,  Evaluation  Ofc,  HDOO 

ATTN:  Library 
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DEPARTMENT  OF  THE  ARMY  (Continued) 

US  Army  Satellite  Comm  Agency 
ATTN:  Doc  Con 

US  Army  TRADOC  Sys  Analysis  Actvy 


ATTN: 

ATAA-PL 

ATTN: 

ATAA-TDC 

ATTN: 

ATAA-TCC,  F.  Payan,  Jr 

US  Army  White  Sands  Missile  Range 

ATTN: 

STEWS-TN-N,  K.  Cummings 

USA  Missile 

Command 

ATTN: 

DRSMI-YSO,  J.  Gamble 

DEPARTMENT  OF  THE  NAVY 

Joint  Cruise  Missiles  Project  Ofc 

ATTN 

JCMG-707 

Naval  Air  Systems  Cmd 

ATTN 

PMA  271 

Naval  Electronic  Systems  Cmd 

ATTN 

PME  106-4,  S.  Kearney 

ATTN 

PME  117-20 

ATTN 

Code  3101,  T.  Hughes 

ATTN 

PME  106-13,  T.  Griffin 

ATTN 

PME  117-2013,  G.  Burnha 

ATTN 

Code  501A 

ATTN 

PME  117-211,  B.  Kruger 

Naval  Intelligence  Support  Ctr 

ATTN 

N I  SC - 50 

Naval  Ocean  Systems  Ctr 

ATTN 

:  Code  532 

ATTN 

Code  5323,  J.  Ferguson 

ATTN 

:  Code  5322,  M.  Paulson 

Naval  Research  Lab 

ATTN 

:  Code  4780 

ATTN 

:  Code  4780,  S.  Ossakow 

ATTN 

:  Code  4187 

ATTN 

:  Code  7500,  8.  Wald 

ATTN 

Code  4700 

ATTN 

:  Code  4720,  J.  Davis 

ATTN 

:  Code  7950,  J.  Goodman 

ATTN 

Code  6700 

Naval  Space  Surveillance  System 

ATTN 

J.  Burton 

Naval  Surface  Weapons  Ctr 
ATTN:  Code  E31 

Naval  Telecommunications  Cmd 
ATTN:  Code  341 

Office  of  the  Deputy  Chief  of  Naval  Ops 
ATTN:  NOP  981 N 

ATTN:  NOP  654,  Strat  Eval  4  Anal  8r 
ATTN:  NOP  9410 

Office  of  Naval  Rsch 

ATTN:  Code  414,  G.  Joiner 
ATTN:  Code  412,  W.  Condell 


DEPARTMENT  OF  THE  NAVY  (Continued) 

Strategic  Systems  Project  Office 
ATTN:  NS  P-2 1 4 1 
ATTN:  NSP-43 
ATTN:  NSP-2722 

Theater  Nuclear  Warfare  Prj  Office 
ATTN:  PM-23,  D.  Smith 

DEPARTMENT  OF  THE  AIR  FORCE 

Air  Force  Geophysics  Lab 

ATTN:  OPR,  H.  Gardiner 
ATTN:  OPR-1 
ATTN:  LKB,  K.  Champion 
ATTN:  CA,  A.  Stair 
ATTN:  PHY,  J.  Buchau 
ATTN:  R.  Babcock 
ATTN:  R.  O'Neil 

Air  Force  Tech  Applications  Ctr 
ATTN:  TN 

Air  Force  Weapons  Lab 
ATTN:  SUL 
ATTN:  NTYC 
ATTN:  NTN 

Air  Force  Wright  Aeronautical  Lab/AAAD 
ATTN:  W.  Hunt 
ATTN:  A.  Johnson 

Air  Logistics  Cmd 

ATTN:  OO-ALC/MM 

Air  University  Library 
ATTN:  AUL-LSE 

Assistant  Chief  of  Staff 

Studies  &  Analyses 

ATTN:  AF/SASC ,  W.  Kraus 
ATTN:  AF/SASC,  C.  Rightmeyer 

Ballistic  Missile  Office/DAA 
ATTN:  ENSN 
ATTN:  SYC,  D.  Kwan 
ATTN:  ENSN,  W.  Wilson 

Deputy  Chief  of  Staff 

Rsch,  Dev  &  Acq 

ATTN:  AFRDP 

ATTN:  AFROS,  Space  Sys  &  C3  Dir 

Deputy  Chief  of  Staff 

Plans  4  Operations 
ATTN:  AFXOKCD 
ATTN:  AFXOKS 
ATTN:  AFXOKT 


Headquarters 
Electronic  Systems  Div 


ATTN: 

ESD/SCTE,  J.  Clark 

ATTN: 

OCT-4, 

J.  Deas 

ATTN: 

SCS-2, 

G.  Vinkels 

ATTN: 

SCS-1E 

68 


DEPARTMENT  OF  THE  AIR  FORCE  (Continued) 


DEPARTMENT  OF  ENERGY  CONTRACTORS 


Foreign  Technology  Div 

ATTN:  TQTD,  B.  Ballard 
ATTN:  Nils,  Library 

Rome  Air  Development  Ctr 

ATTN:  OCS,  V.  Coyne 
ATTN:  TSLD 

Rome  Air  Development  Ctr 

ATTN:  EEP,  J.  Rasmussen 

Headquarters 
Space  Command 

ATTN:  DC,  T.  Long 

Space  Div 

ATTN:  YGJB,  W.  Mercer 
ATTN:  YKM ,  E.  Norton 
ATTN:  YKM,  W.  Alexander 

Strategic  Air  Command 
ATTN:  NRT 

ATTN:  OCXT,  T.  Jorgensen 
ATTN:  XPFS 
ATTN;  DCS 

ATTN:  ADWAT,  R.  Wittier 

DEPARTMENT  OF  ENERGY 

Department  of  Energy 
ATTN:  DP-233 

OTHER  GOVERNMENT  AGENCIES 

Central  Intelligence  Agency 
ATTN:  OSWR/NED 

ATTN:  OSWR/SSD  for  K.  Feuerpfetl 

Department  of  Commerce 
National  Bureau  of  Standards 

ATTN;  Sec  Ofc  for  R.  Moore 

Department  of  Commerce 
National  Oceanic  &  Atmospheric  Admin 
ATTN:  R.  Grutb 

National  Telecommunications  Sciences 
ATTN:  L.  Berry 
ATTN:  A.  Jean 
ATTN:  W.  Utlaut 

Department  of  State 

Office  of  International  Security  Policy 
ATTN:  PM/ STM 

NATO 

NATO  School,  SHAPE 

ATTN:  US  Documents  Officer 


EG&G,  Inc 

ATTN:  J.  Colvin 
ATTN:  D.  Wright 

University  of  California 

Lawrence  Livermore  National  Lab 
ATTN:  L-389,  R.  Ott 
ATTN:  L-31,  R.  Hager 
ATTN:  Tech  Info  Dept,  Library 

Los  Alamos  National  Lab 

ATTN:  MS  670,  J.  Hopkins 
ATTN:  P.  Keaton 
ATTN:  MS  664,  J.  Zinn 
ATTN:  T.  Kunkle,  ESS-5 
ATTN:  R.  Jeffries 
ATTN:  D.  Simons 
ATTN:  J.  Wolcott 

Sandia  National  Labs 

ATTN:  B.  Murphey 
ATTN:  T.  Cook 

Sandia  National  Labs 

ATTN:  D.  Thornbrough 
ATTN:  Tech  Library,  3141 
ATTN:  D.  Dahlgren 
ATTN:  Space  Project  Div 
ATTN:  Org  4231,  T.  Wright 
ATTN:  Org  1250,  W.  Brown 

DEPARTMENT  OF  DEFENSE  CONTRACTORS 

Aerospace  Corp 

ATTN:  V.  Josephson 
ATTN:  R.  Slaughter 
ATTN:  T.  Salmi 
ATTN:  I.  Garfunkel 
ATTN:  K,  Cho 
ATTN:  J.  Straus 
ATTN:  D.  Olsen 

Analytical  Systems  Engrg  Corp 
ATTN:  Radio  Sciences 

Analytical  Systems  Engrg  Corp 
ATTN:  Security 

BDM  Corp 

ATTN:  T.  Neighbors 
ATTN:  L.  Jacobs 

Berkeley  Rsch  Associates,  Inc 
ATTN:  C.  Prettie 
ATTN:  J.  Workman 
ATTN:  S.  Brecht 

Boeing  Aerospace  Co 

ATTN:  MS/87-63,  D.  Clauson 
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DEPARTMENT  OF  DEFENSE  CONTRACTORS  (Continued) 

DEPARTMENT 

OF  DEFENSE  CONTRACTORS 

Boeing  Co 

HSS,  Inc 

ATTN:  G.  Hall 

ATTN:  S.  Tashird 

ATTN 

D.  Hansen 

IBM  Corp 

BR  Communications 

ATTN:  J.  McLaughlin 

ATTN 

H.  Ulander 

Institute  for  Defense  Analyses 

University  of  California  at  San  Diego 

ATTN 

J.  Aein 

ATTN:  H.  Booker 

ATTN 

E.  Bauer 

ATTN 

F.  Wolfhard 

Charles  Stark  Draper  Lab,  [nc 

ATTN 

H.  Gates 

ATTN:  D.  Cox 

ATTN:  A.  Tetewski 

International  Tel  Telegraph  Corp 

ATTN:  J.  Gilmore 

ATTN 

Technical  Library 

Computer  Sciences  Corp 

International  Tel  &  Telegraph  Corp 

ATTN:  F.  Eisenbarth 

ATTN 

G,  Wetmore 

Comsat  Labs 

JAVCOR 

ATTN:  G.  Hyde 

ATTN:  D.  Fang 

ATTN 

J.  Sperling 

Johns  Hopkins  University 

Cornell  University 

ATTN 

T.  Evans 

ATTN:  D.  Farley,  Jr 

ATTN 

K.  Potocki 

ATTN:  M.  Kelly 

ATTN 

J.  Newland 

ATTN 

P.  Komiske 

E-Systems,  Inc 

ATTN 

J.  Phillips 

ATTN:  R.  Berezdivin 

Kaman  Tempo 

Electrospace  Systems,  Inc 

ATTN 

B.  Gambill 

ATTN:  H.  Logston 

ATTN 

K.  Schwartz 

ATTN:  P.  Phillips 

ATTN 

DAI  SAC 

ATTN 

J.  Devore 

EOS  Technologies,  Inc 

ATTN 

W.  McNamara 

ATTN:  8.  Gabbard 

Kaman  Tempo 

General  Electric  Co 

ATTN 

DAS  I  AC 

ATTN:  A.  Steinmayer 

ATTN:  C.  Zierdt 

Litton  Systems,  Inc 

ATTN 

B.  Zimmer 

General  Electric  Co 

ATTN:  G.  Millman 

Lockheed  Missiles  &  Space  Co,  Inc 

ATTN 

J .  Kumer 

General  Rsch  Corp 

ATTN 

R.  Sears 

ATTN:  B.  Bennett 

Lockheed  Missiles  S  Space  Co,  Inc 

GEO  Centers,  Inc 

ATTN 

Dept  60-1Z 

ATTN:  E.  Marram 

2  cy  ATTN 

D.  Churchil 1 ,  Dept  62 

GTE  Communications  Products  Corp 

M.I.T.  Lincoln  Lab 

ATTN:  R.  Steinhoff 

ATTN 

D.  Towle 

GTE  Communications  Products  Corp 

MA/COM  Linkabit  Inc 

ATTN:  J.  Concordia 

ATTN 

A.  Viterbi 

ATTN:  I.  Kohl  berg 

ATTN 

H.  Van  Trees 

ATTN 

I .  Jacobs 

Harris  Corp 

ATTN:  E.  Knick 

Magnavox  Go^t  &  Indus  Electronics 

ATTN 

G.  White 

Honeywell,  Inc 

ATTN:  G.  Terry,  Avionics  Dept 

McDonnell  Douglas  Corp 

ATTN:  A.  Kearns,  MS924-3 

ATTN 

Tech  Library  Svcs 

ATTN 

W.  Olson 

Horizons  Technology,  Inc 

ATTN 

R.  Halprin 

ATTN:  R.  Kruger 
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DEPARTMENT  OF  DEFENSE  CONTRACTORS 

Meteor  Communications  Corp 

Rand  Corp 

ATTN 

R.  Leader 

ATTN 

E.  Bedrozian 

ATTN 

C.  Crain 

Mission  Rsch  Corp 

ATTN 

P.  Davis 

ATTN 

C.  Lauer 

ATTN 

Tech  Library 

Riverside  Rsch  Institute 

ATTN 

F.  Fajen 

ATTN 

V.  Trapani 

ATTN 

R.  Kilb 

ATTN 

S.  Gutsche 

Rockwell  International  Corp 

ATTN 

R.  Bigoni 

ATTN 

R.  Buckner 

ATTN 

R.  Bogusch 

ATTN 

F.  Guigliano 

Rockwell  International  Corp 

ATTN 

G.  McCartor 

ATTN 

S.  Quil ici 

ATTN 

R.  Hendrick 

4  cy  ATTN 

0.  Knepp 

Science  Applications,  Inc 

5  cy  ATTN 

Doc  Con 

ATTN 

L.  Linson 

ATTN 

D.  Hamlin 

Mitre  Corp 

ATTN 

C .  Smith 

ATTN 

MS  J104,  M.  Dresp 

ATTN 

E.  Straker 

ATTN 

A.  Kymmel 

ATTN 

G.  Harding 

Science  Appl ications,  Inc 

ATTN 

C.  Callahan 

ATTN 

J.  Cockayne 

Mitre  Corp 

Science  Appl ications,  Inc 

ATTN 

W.  Hall 

ATTN 

M.  Cross 

ATTN 

J.  Wheeler 

ATTN 

M.  Horrocks 

SRI  International 

ATTN 

W.  Foster 

ATTN 

R.  Tsunoda 

ATTN 

J.  Vickrey 

Pacific-Sierra  Rsch  Corp 

ATTN 

W.  Chesnut 

ATTN 

E.  Field,  Jr 

ATTN 

R.  Leadabrand 

ATTN 

F.  Thomas 

ATTN 

R.  Livingston 

ATTN 

H.  Brode,  Chairman  SAGE 

ATTN 

D.  McDaniels 

ATTN 

M.  Baron 

Pennsylvania  State  University 

ATTN 

G.  Price 

ATTN 

Ionospheric  Research  Lab 

ATTN 

D.  Neil  son 

ATTN 

A.  Burns 

Photometries,  Inc 

ATTN 

W.  Jaye 

ATTN:  I.  Kofsky 

ATTN 

J.  Petrickes 

ATTN 

C.  Rino 

Physical  Dynamics,  Inc 

ATTN 

V.  Gonzales 

ATTN 

E.  Fremouw 

ATTN 

G.  Smith 

Physical  Rsch,  Inc 

Stewart  Radiance  Lab 

ATTN 

R.  Deliberis 

ATTN:  J.  Ulwick 

ATTN:  T.  Stephens 

Technology  International  Corp 

RAD  Associates 

ATTN 

W.  Boquist 

ATTN 

R.  Turco 

ATTN 

W.  Karzas 

Toyon  Rsch  Corp 

ATTN 

H.  Ory 

ATTN 

J.  Ise 

ATTN 

C.  Grei finger 

ATTN 

J.  Garbarino 

ATTN 

P.  Haas 

ATTN 

M.  Gantsweg 

Tri-Com,  Inc 

ATTN 

F.  Gilmore 

ATTN 

D.  Murray 

ATTN 

W.  Wright 

TRW  Electronics  &  Defense  Sector 

R&D  Associates 

ATTN 

R.  Plebuch 

ATTN 

B.  Yoon 

ATTN 

G.  ICirchner 

Rand  Corp 

ATTN:  B.  Bennett 
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Utah  State  University 
Attention  Sec  Con  Ofc  for 


ATTN: 

K.  Baker, 

Dir  Atmos  &  Space  Sci 

ATTN: 

L.  Jensen 

Elec  Eng  Dept 

ATTN: 

D.  Burt 

ATTN: 

A.  Steed 

DEPARTMENT  OF  DEFENSE  CONTRACTORS  (Continued) 
VisiDyne,  Inc 


ATTN: 

C.  Humphrey 

ATTN: 

0.  Shepard 

ATTN: 

W.  Reidy 

ATTN: 

J.  Carpenter 
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